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I study the incorporation of renormalization group (RG) improved BFKL kernels in the 
' Balitsky-Kovchegov (BK) equation which describes parton saturation. The RG improve- 

ment takes into account important parts of the next-to-leading and higher order logarithmic 
corrections to the kernel. The traveling wave front method for analyzing the BK equation 
| is generalized to deal with RG-resummed kernels, restricting to the interesting case of fixed 

O ■ QCD coupling. The results show that the higher order corrections suppress the rapid in- 

crease of the saturation scale with increasing rapidity. I also perform a "diffusive" differential 
(-H ■ equation approximation, which illustrates that some important qualitative properties of the 

| kernel change when including RG corrections. 

> : 

m . 

O ■ I. INTRODUCTION 

O ' 

<n : 

y—( \ Studies of parton saturation and the unitarization of high energy amplitudes began more than 

20 years ago with the work of Gribov, Levin and Ryskin [l| and has now evolved into a very 
active field of research Q]. The simplest modern QCD evolution equation that describes parton 
saturation in the high energy limit is the Balitsky-Kovchegov (BK) equation 0, Q]. The BK 
equation is a "mean-field" approximation which neglects higher order correlators present in the so- 
called Balitsky [4] and JIMWLK equations (see j2] and references therein). Nevertheless, because 
the corrections neglected in the mean field approximation are likely small at phenomenologically 
interesting energies jH], it is expected to be accurate for dipole scattering on large nuclei and at 
least a phenomenological model for saturation in scattering on protons. The equation has the 
correct high-energy limit at small parton densities, when it reduces to the BFKL equation [jj], but 
it has a non-linear term which reduces the power-like growth with energy. 

The BK equation is derived at leading logarithmic (LL) accuracy and thus reduces to the LL 
BFKL equation in the linear regime. The next-to-leading logarithmic (NLL) corrections to the 
BFKL equation were computed in 0j and were found to be negative and very large. As the BFKL 
kernel is a central part of the BK equation, it is important to understand the effect of the NLL 
corrections on the solutions of the BK equation. 

In fact the NLL BFKL kernel is larger than the leading order kernel, and for any reasonable 
value of a s it is negative; the pomeron intercept becomes less than one for a s > 0.15. The 
kernel also has two complex conjugate saddle points which lead to oscillating cross sections. The 
perturbative expansion of the BFKL kernel is therefore highly unstable and far from converging. 
These problems are discussed in detail in [8|]. 

A cure was proposed by Ciafaloni, Colferai, and Salam [{J, who realized that the large corrections 
come from collinear contributions of the DGLAP type. They devised the so-called w-resummation, 
which makes the NLL kernel compatible with renormalization group (RG) requirements through 
matching to the DGLAP limit and resummation of spurious poles. A somewhat different scheme 
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has been proposed in Many groups have also studied ways of including some of the most 

important NLL effects in the LL kernel [0, E3, A common feature of these schemes is that 
they modify the functional dependence of the Mellin space BFKL kernel x(l) so that it acquires 
a dependence also on the to variable. This will be discussed at length below. In a slight abuse of 
terminology I will henceforth refer to all these schemes as "RG improved" kernels. 

Another important part of the NLL corrections is the running of the QCD coupling. It is, 
however, not clear how to incorporate the effect of the running coupling in the BK and Balitsky- 
JIMWLK equations, and recent studies have come to different conclusions LJ]. It is also reasonable, 
but by no means obvious, to assume that the NLL corrections to the BFKL part of the BK equation 
will be the same as for the BFKL equation. Since the way to include the running coupling is unclear 
it is therefore useful to study the BK equation with NLL corrections to the BFKL kernel but with 
a fixed coupling. 

The onset of saturation corrections is governed by the saturation scale Q s — the momentum that 
separates the regions dominated by linear evolution and saturation. Because of the non-linear term 
in the BK equation it is difficult to obtain analytical solutions, but at large momenta the amplitude 
is small enough that one may consider the linear equation only (the BFKL equation) together with 
appropriate boundary conditions. This has been used in a series of papers 0, EllzI, El] to 
compute the rapidity dependence of Q s = Q S (Y). Recently the method of traveling waves has 
been developed 19, 2(], 21], which allows performing these calculations in a more systematic way. 
Asymptotically the saturation scale increases exponentially as rapidity increases 

EH ED, 

but, as 

we will discuss below, there is also a factor with a sub-asymptotic Y-dependence which slows down 
the increase 13, LL8|, |l9( . 

Saturation is intimately related to geometric scaling, which was discovered in the HERA data 
in Ref. [22I]. It refers to a scaling of the j*p cross section for small Bjorken-x, where the cross 
section cr 7 p (x,Q 2 ) is not a function of x and Q 2 separately, but only through the combination 
r = Q 2 /Q 2 s (x). As we shall see, the BK equation predicts geometric scaling of the amplitude and 
a specific pattern of violation of this scaling. 

In this paper I make an analytical study using the method of traveling wave fronts to compute 
the higher-order corrections and suggest a method which is simple to implement for any higher 
order corrected kernel (see also 23f] ) . I will do this for the case of a fixed strong coupling a. The 
case of an RG-improved kernel with running coupling was very recently done by Peschanski and 
Sapeta [24(] and will therefore not be studied here. The present paper should therefore be taken 
together with 2J] for a treatment of the RG-improved BK equation. 

The outline of the paper is as follows. Section [II] contains a brief review of the traveling wave 
front method for analyzing the BK equation and section IIIII discusses higher order corrections to 
the BFKL kernel. In Section IIVI I then generalize the method described in Section [TT] to deal with 
RG-resummed kernels. I show how to compute analytic results for the saturation scale and the 
shape of the solution in momentum space above the saturation scale, restricting to the case of fixed 
coupling. This is developed for a general form of the resummed kernel and then applied to various 
schemes proposed in the literature in Section [Vj In Section I VII finally, I perform a "diffusive" 
differential equation approximation, which illustrates that some important qualitative properties 
of the kernel change when including RG corrections. 



II. THE BK EQUATION AND TRAVELING WAVES 



The BK equation in position space describes the scattering amplitude A/"(xoi, Y) of a dipole on 
a target (originally a large nucleus), where the target is specified by the initial condition. For large 
and homogeneous targets (target radius much larger than the involved dipoles in the cascade) one 
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may neglect the dependence on the impact parameter and Fourier transform to momentum space. 
The equation can then be written on the form [25| 



dyM = a X o(-d L )M - aN 2 . 



(1) 



where at = N c a s /ir, L = h\{k 2 /k^), with ko some low momentum scale, e.g., Aqcd, and dy = 
d/dY, di = d/dL. The quantity jV(k, Y) which we will study in the following is related to the 
unintegrated gluon distribution of the target. The non-linear term in Eq. ([I]) is particularly simple 
in this representation and is not acted on by the kernel operator. Linearizing Eq. ([1]) one arrives 
at the BFKL equation [g]. 

The operator xo(~ 9l) is given by the characteristic BFKL function 

X o(i) =2^(l)-V>(7)-V>(l-7), (2) 

the Mellin transform of the LL BFKL kernel in momentum space 0]. The explicit form of the 
term xo(~ is a convolution Kbfkl <8> of the BFKL integral kernel -Kbfkl and the gluon 
distribution, but it can be formally represented as an infinite order differential operator through 
the power series 



Xo(-d L )M: 



Xo(7o) + {-0 L - 7o)x'o(7o) + 2\(.-°L ~ lofx'oilo) + • 



(3) 



where 70 is a constant. 

The approximations involved in deriving the BK equation are (i) large iV c , (ii) no correlations 
between different dipoles, and (hi) leading logarithmic (LL) resummation. In this paper we will 
investigate the corrections to the third point. The second point has been studied a lot in the last 
years, e.g., in litSj]. 

Munier and Peschanski 19, 2(], 21] found an interesting relationship between non- linear evolu- 
tion equations in QCD (in particular the BK equation) and statistical physics studies of traveling 
wave fronts. These are wave fronts that separate a stable state from an unstable state, which 
invades the unstable state by propagating in space with an essentially fixed shape. Such fronts are 
described by non-linear diffusion equations with a growth term and a reduction term 27, 28|, 29, 30]. 

They observed [19|] that if one truncates the series ([3|) and keeps only the terms up to order d\ 
one obtains the differential equation 



dyM 



a 



Xo(7o) + (-01 



7oko(7o) + 2\(- d L 



7o) 2 Xo(7o) 



M - 6tN 2 



(4) 



which, with a suitable change of variables Y oc t and L 
gives the equation 



Ax + Bt (with A,B some constants), 



d t u 



d x u + u 



(5) 



This is the celebrated Fisher-Kolmogorov-Petrovsky-Piscunov (FKPP) equation 27|, |28|, which 
under certain conditions leads to traveling wave solutions u = u(x,t). The condition is that 
the initial conditions are steep enough, i.e. decay like exp(— Ax) at large x with A large enough, 
see [2^, [13] for more details. I will not discuss the initial conditions further here, but will assume 
that they are steep enough, as is the case in QCD [191 ]. 

The traveling wave property means essentially that the solutions are saturated wave fronts in 
the space variable x, which propagate forward in space as the time t increases, with a fixed (or 
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nearly fixed) shape. The speed of propagation v(t) of this front has been known for a long time 311 ] : 
it approaches an asymptotic speed v* like 

«(*) = - ^ + ^r 3/2 ) (6) 



where for the FKPP equation v* = 2 and A* = 1. The result for the saturation scale given in [It 
is obtained from this formula by using the correct change of variables to get the approximate BK 
equation in the form of the FKPP equation, Eq. ([5]). That is, when studying the BK equation 
in the traveling wave framework, we study the evolution of a front that propagates in the "space" 
variable L ~ lnfc 2 as the "time" Y increases. 

The 0(t~ 3 / 2 ) term in the velocity relation for FKPP-like equations has been found by Ebert 
and van Saarloos [30], namely 



where A* is related to the wave number and has the dimensions of inverse length, and D is the 
diffusion constant. (For the FKPP equation ([5]), D = 1.) It has also been shown that these are all 
universal terms in the velocity relation, i.e., all other terms depend on the initial conditions. The 
three constants v*, A* and D completely determine the velocity relation of the traveling waves, and 
thus the saturation scale in the case of the BK equation. These constants are calculated from the 
dispersion relation for Fourier modes of the solutions of the equation under study. 

The most important point of the approach of Munier and Peschanski may be that the velocity 
that one calculates is independent of the exact form of the non-linearity, i.e., it is enough to analyze 
the linearized integro-differential equation, i.e., the BFKL equation, which of course is much easier 
than the full BK equation. This does not mean that one gets saturation from the linear equation — 
the non-linearity is crucial in selecting the velocity (J7|). See Section III of Ref. [30] for a thorough 
discussion of this point. 



III. HIGHER ORDER CORRECTIONS TO THE BFKL AND BK EQUATIONS 

The full NLL BFKL kernel is given by x{l) = Xo(7) + c^Xiil) where xo is the LL BFKL kernel 
and xi is the second order contribution to the BFKL kernel and is given in Refs. 0. The full 
BK equation has not been derived at NLL accuracy — only the BFKL equation is known. There 
will also be corrections to the non-linear term, but it is reasonable to assume that the linearized 
equation will be the same as the NLL BFKL equation. Since the observables studied here are, 
as explained above, determined by the linear equation we can study the BK equation with NLL 
corrections only in the linear term and obtain the correct Independence of the saturation scale. 
The absolute normalization of Q s is not constrained by the method. The normalization of the wave 
front, however, is given by the full non- linear evolution of the initial condition and may change if 
the non-linear term is changed. The same holds for the wave front in the saturation region, but 
here we only aspire to calculating it in the region above the saturation scale. 

Using the improved kernels mentioned in the Introduction, the BFKL kernel is modified to the 
form x(7, to) depending on both 7 and to. The RG-improved BK equation is then non-local in both 
L and Y and mixes transverse and longitudinal degrees of freedom (see also Section IVI|) . It can be 
written 



8yN = a X {-d L , c\)N - aN 2 



(8) 
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regardless of the specific form of the kernel. Here one can take either fixed or running coupling. 
The evolution equation is now non-local in rapidity, whereas the uncorrected kernel (LL or NLL) 
only is non-local in momentum. This is where the traveling wave method reveals its power — it can 
easily be generalized to kernels much more complicated than LL BFKL (and similarly the FKPP 
equation can be turned into higher order differential equations). 

The first study of hi gher order corrections to the saturation scale was to my knowledge done 
by Golec-Biernat et al. [32j in their numerical study of the BK equation, by implementing the 
kinematical constraint of Ref. The first analytical study was performed by Triantafyllopoulos 
[l~8| using BFKL evolution with absorbing boundary conditions [13]. The higher corrections in 
this case were taken into account through the use of the resummed NLL kernel of Refs. 0] with a 
running coupling constant. Khoze et al. [l3| have made another BFKL study where they compute 
the asymptotic term in the saturation scale using a resummed LL kernel along the lines of i.e., 
a simpler variant than the one in 18[] . Chachamis et al. 33f] used the BFKL kernel with a rapidity 
veto [12] in the BK equation and performed a numerical solution of this BK equation. Finally, 
Gotsman et al. 34] studied a DGLAP-corrected BFKL kernel similar to the one used in 0. All 
these studies naturally find a reduced growth with energy of the saturation momentum. 



IV. METHOD FOR CALCULATION 
A. The leading logarithmic BK equation 

The linearized equation corresponding to (pQ) is the LL BFKL equation, which can be solved 
analytically by Mellin transformation; the solution is given by the double inverse Mellin transform 1 

m,Y) = / fr / ^Mo{lM- — ry (9) 

where Mn[j,uj) encodes the initial conditions. This is analogous to the Fourier-Laplace transform 
used in [30J to analyze the FKPP equation leading to a dispersion relation between wave number 
and frequency for the Fourier modes 2 . In the present case we obtain a dispersion relation for Mellin 
moments by picking up the pole in the denominator to perform the w-integral. This gives 

M(k,Y) = [ ^lA/ (7)e-( 1 -^ L+ ^°W y . (10) 
Jc 27ri 

We thus have the dispersion relation u)(j) = axo(l) for the solution of the BK equation. The 
evolution variable is now rapidity Y, which corresponds to time t in the FKPP equation. 

Let us now see how easily the saturation scale is computed using the Ebert-van Saarloos (EvS) 
method [3p| . The calculation consists of finding the three parameters v* , A* = k* = 1 — j c , and D. 
The linearized LL-BK equation (i.e., the LL-BFKL equation) is in Mellin 7-space written as 

dyfa,Y) = a X0 (7)0(7, Y) (11) 

where 0(7, Y) denotes the Mellin transform of M(k, Y) with respect to L. Following the general 
formalism of EvS Section V-D, we write this as 

d Y + J( 1 )]^( 1 ,Y) = (12) 



1 Note that I use a slightly different convention than e.g. [19( ] for the 7-dependence of the exponent. 

2 Working with Mellin transforms instead of Fourier-Laplace transforms, the quantities uj and k* are real instead of 
purely imaginary. Appropriate changes must therefore be made to the formulas of [30| . 
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where — Tcj) is the Mellin transform with respect to k 2 of -Kbfkl <8> M, i.e., T = —0x0(7)- Define 
also S(7,a;), the Mellin transform with respect to Y of the left hand side of (|12p . i.e., it is the 
double Mellin transform of the BFKL equation. In general, for differential equations, S and T are 
M x M matrices where M is the order of the highest time derivative. 

The dispersion relation for Mellin moments is obtained from the requirement £(7, u) = where 
S is the characteristic function 5(7, u) = det §(7, u) = 0x0(7)- This again yields 0^(7) = «Xo(7) 
as expected. One now finds the asymptotic velocity v* and the saddle point k* = X* = 1 — 7 C by 
solving the set of equations [EvS Eqs. (5.60) and (5.66)] 

<Mt) 



w(7c) 



7 - ( 



(13) 



1 - 7c 1 - 7c 

The relevant saddle point for evolution with constant amplitude is therefore not the usual BFKL 
saddle point 70 = 1/2 but is given by the solution of the equation Xo(7c)(l — 7c) + Xo(7c) = 
as was shown already in A numerical solution of (|13|) gives 7 C = 0.373. Finally, the diffusion 
coefficient D is obtained from 

1 d 2 uo{~i) 



D 



2! dj 2 



a 



Xo(7c)- 



7c 



The speed of the front is now given by Eq. ([7]) and the position is then simply 



x(Y) = \n[Q 2 (Y)/k 2 } = v*Y - JL In Y - ^JL 



7T 

DY 



which, inserting v* , A*and D gives the result (neglecting ko) 



lnQ 2 (Y) 



a- 



. Xo(7c 
l-7o 



■Y 



2(1 - 7c 



hxY 



2tt 



(l-7c) 2 Y«Xo / (7c)^ 



+ 0{l/Y). 



(14) 



(15) 



(16) 



This expression for the saturation scale of a system whose evolution is described by Eq. (pQ) was 
first obtained in 2l|] . These three terms are universal. The first term, which gives the asymptotic 
exponential rise of the saturation scale used e.g. in the saturation model of [3511 . was obtained 
already in [l] from the double logarithmic approximation, and from LL BFKL in jig ]. The second 
term was obtained in 17| using BFKL evolution with an absorbing boundary to mimic the non- 
linearities and was confirmed in 1{^, 20] using the traveling wave method. 3 

The shape of the wave front in the so-called lea ding edge region can also be computed in the 
traveling wave framework. Ebert and van Saarloos [30| compute the solution at the same level of 
approximation as the velocity ([7]), i.e., including all universal terms. Brunet and Derrida [36| use 
a somewhat simpler method which allows computing the solution at the level of approximation of 
([6]), which gives simpler results. This method was employed in Ref. [2ol ] to compute the solution of 
the LL BK equation in the leading edge region for both fixed and running coupling. The result is 

—1 



M{k,Y) 



In 



QtOOJ \Q 2 (Y) 



exp 



1 



ADY 



■In 1 



In 



k 2 



«X '(7c) \Q 2 (Y)J \Q 2 (Y) 



exp 



Qs( Y ) 
1 



(17) 



hr 



2«x"(7c)^ \Q 2 (Y) 



3 Note that when comparing the results of pjj and [l9l . I20I ] there is a difference of a factor of 2 due to the convention 
for the definition of the BFKL characteristic function xo- Furthermore, there are some terms in the results of [13] 
that are non-universal and subleading in the classification used here. 
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up to a constant factor. This result has been obtained both using an absorbing boundary [17]] and 
traveling waves (2o| . 

The agreement between these completely different methods both for the saturation scale and 
the wave front shape is good evidence that the travelingwave method is reliable and general. 
Furthermore, the numerical study of the BK equation in [5[ shows very good agreement at large 
Y with Eq. (|16p and the corresponding analytical formulas for J\T(k, Y) in the leading edge region. 

The description in terms of operators S and T might seem unnecessarily formal, but a main 
result of the analysis of EvS [30| is that the method generalizes to much more general types of 
equations than differential equations of the FKPP type, 4 for example higher order differential 
equations, integral equations with a kernel in either time (memory kernels) or in space (non-local 
kernels), or even difference equations. 

Thus, the kernel in momentum or configuration space might have any complicated structure, 
but if one knows the Mellin (or Fourier) transform of the kernel, one can solve for the dispersion 
relation and compute v*, A* and D (in fact, in the following treatment the higher order corrections 
are implemented in Mellin space rather than in momentum space). This dispersion relation can 
lead to non-algebraic equations, however, unlike for the FKPP equation. This generalization will 
be exploited below for the RG-resummed case. 



B. RG-improved BK 



Let me now generalize to RG-improved kernels at fixed coupling. The dispersion relation now 
reads 



u(7,a) = ax(7>^(7,a)) 



(18) 



which implicitly defines the function lo. For each value of 7 and a Eq. (|18p gives a transcendental 
equation, which can in principle be solved numerically for 0^(7). This is not what we will do here. 
The asymptotic velocity v* and saddle point A* = 1 — 7 C are as earlier found from the equations 



do; (7, a) 



l-7c ' 



7c 



(19) 
(20) 



but since cj(j,a) does not have an explicit form it is not as straightforward to obtain the required 
quantities. For LL there is a saddle point of the characteristic function Xo(7) determined from the 
equation given above, but now the function x = x{li w ) nas an additional variable. This means 
that we will have to determine two constants ~f c and iv c . To keep notation simple, let us introduce 
the shorthand Xc = x(lc,u c ), x' c = djXilc^c) and Xc = d^xilc^c)- 

First, we require uj c = ax(7c,w c ). Observing that on one hand, v* is given by 



did 

d'y 



-a 



7c 



dx(7,u) dx(j,u) duj 



97 



duj c?7 



7c 



so that 



dj 



ax'c 
1 - atXc 



(21) 



(22) 



4 They have to be equivalent in the sense that they must be reaction-diffusion equation with second-order diffusion. 
The running coupling case, for example, leads to a different equivalence class which has D — 0. 
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and using Eq. ([20]) we obtain — x£(l — 7 C ) = Xc [1 — «Xc] = Xc — w c Xc- We therefore have the 
system of equations 

jxc(! -7c) +Xc- UcXc = 
\w c = «Xc 

The first equation reduces to the corresponding LL one if there is no dependence on u in the kernel. 

The advantage of knowing both 7 C and uj c is obvious — there exists a closed form for all the 
constants that will appear in the expression for the saturation scale, and it is not necessary to use 
functions that can only be defined numerically and to numerically differentiate these functions. It 
is much more straightforward to just numerically solve (I23h . 

Once we have determined 7 C and u c , the other constants v* and D are fixed (remember A* = 
1 — 7c)- We are then in a position to write down the expression for the velocity of the front, 



w-'-w^V-y/trwsF) (24) 

where v* = u> c /(l — j c ) and D is obtained by applying the chain rule as in Eq. ([2~T1) twice; 



D = - 



i J2 / / \ 2 



2 d-y 



7c 



2x'c(7c-l) r c l-7c Vl-7 



x' c ' + 2-^-x'c+ 7^ Xc • (25) 



Finally this gives the saturation scale, 



2 



-1/2 



m Ql(Y) = ^-Y - In r - J 18 ^ 3y U + iJU^Z+fjS!-) ic] 

l-7c 2(1 -7c) \j u c {-f c - l^Y y l- 7c \l-7cj J 

(26) 

All the discussion so far is perfectly general in that we have not chosen any particular form for 
the NLL BFKL kernel, we have only anticipated that it will depend on both 7 and u. 



APPLICATION TO RG-IMPROVED KERNELS 



The method from the last section can now be easily applied to some of the RG-improved kernels 
suggested in the literature. I will look at the rapidity veto model of 12|], the kinematical constraint 



model of and the resummed model of 13|]. I will thus solve Eqs. (|23|) for these models 



numerically using Mathematica, and use the obtained values for 7 c ,a; c in the expression (J26J) for 
the saturation scale and in the generalized expression (|40p below for the shape of the wave front. 



A. Rapidity veto 



The idea of a rapidity veto in the kernel 12|] is that by forbidding gluon emissions that are very 
close in rapidity one accounts for a part of the NLL corrections. This was studied numerically in 
the context of the BK equation in Ref. 



33] 



The form of the kernel is 



X(7,^) = Xo(7)e VU} 



(27) 



where r/ is the minimum rapidity interval imposed on two subsequent emissions, and is taken as a 
parameter. A numerical solution for 0^(7) is shown in Fig. [TJ 
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The equations (J23J) thus take the form 

w c = axo(7c)e-'"" c 



Xohchce-^ = Xo(7 c )e-^ c (l + ar, X o(lc)e~^) 



(28) 



which in fact for this factorized form of the kernel allows eliminating u c from the result so that 
one arrives at the following equation for 7 C : 



x'o(7c)(l-7c) , - — / \ f, , Xo(7c)(l-7c) 
- + 1 = -arixoilc) exp 1 + 



xodc 



xodc 



(29) 



which, taking into account the difference in convention, 7 — > 1 — 7, is the result for 7 C found in 
Ref. [33J]. Note that for r\ — ► this reduces to the LL result. 

The authors of [33J give an analytic result for the saturation scale with only the asymptotic 
term for the velocity and then solve the BK equation numerically. But using the present method 
it is easy to find the subleading terms for the saturation scale. To do this we need to find the 
diffusion coefficient given by Eq. (|25p . 



D 



2Xo(7c)(7 c " 1) 



Xo(7c)-2r,Xo(7c) 



«Xo 7c )e ^ 2 

: + V Xc 7c 

1 -7c 



«Xo(7c)e-^ 
l-7c 



with 



exp 1 + 



Xo(7c)(l ~7c) 
Xo(7c) 



Thus, the saturation scale with all universal terms reads 



m^On = a^Y exp ( 1 + X ' M{1 " 7c 



1 -7c" V Xo(7 
with 7 C and D given by Eqs. (|29p and (|30p . 



2(1 -7c) 



(1 - 7c ) 2 V DY 



7T 



(32) 



B. Kinematical constraint 



The so-called kinematical constraint, or consistency constraint ll|, imposes a cut-off on the 
virtualities of gluons in the real emission part of the BFKL kernel, motivated by the requirement 
that in the high energy kinematics considered, the virtualities of the exchanged gluons should be 
dominated by their transverse momenta. It is also a part of the resummed framework of [9] . 

This constraint has been shown to reproduce a large part of the NLL-corrections to the pomeron 
intercept. It has been used in several phenomenological applications of BFKL dynamics. Because 
the value of the pomeron intercept is reduced, the results show a reduced growth of cross sections 
with energy, and generally fit data better than LL BFKL. 

The constraint implies shifting the 7 in the kernel by ±cj/2, 



X (7, w) = 2^(1) - V (7 + |) " V> (l " 7 + |) 



(33) 



such that, in contrast to the kernel f)2T[) with rapidity veto discussed above, the kernel (j33j) does 
not have poles in 7 = 0, 1. 

In Fig. [T]the solutions of the implicit equation for a; (7, a) for a fixed at = 0.2 for the kernel (|33p 
is compared with the LL kernel xo- I also show the kernel (|27|) with the rapidity veto for a choice 
of r\ that gives the same pomeron intercept as the kinematical constraint. 
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FIG. 1: Characteristic functions 07(7, a = 0.2). The solid line is LL BFKL, Eq. ([2|), the dashed line is the 
kinematical constraint, Eq. ([33]) . and the dotted line is the rapidity veto, Eq. (j27]) . for 77 = 1.87, which gives 
the same value for the pomeron intercept as the dashed line. 



The kernel f|33|) above is symmetric with respect to the interchange 7 — ► 1 — 7, just as the 
original LL BFKL kernel. It is thus written for a symmetric scale choice so in the evolution variable 
Y = m(s/so), were the scales on both sides of the ladder appear in so- In the present analysis we 
would like to make an asymmetric scale choice. In such a case the kinematical constraint should 
be formulated slightly differently [9], namely one writes the kernel as 

X ( 7 ,u;) = 2^(l)-V(7)-V'(l-7 + ^), (34) 
which is now asymmetric with respect to 7 —* 1 — 7. This form will be used below. 



C. Khoze Martin Ryskin Stirling model 

Khoze, Martin, Ryskin and Stirling (KMRS) [131] proposed a simpler model for the NLL cor- 
rected LL kernel, which is based on the approach of [9( but with a fixed coupling. This kernel has 
also been used to estimate the effect of higher order corrections in diffractive 7*7* — > pp production 



37J which provides a possibility to test it against the full NLL calculation [381 ] . 
The characteristic function is expressed as 

X(7>^) = Xo(7) + «Xi(7>^) (35) 

where Xo{l) is the usual BFKL function ([2]). The correction piece Xi{li w ) is given by 

_ 1 + aMiM 1 1+faMiM 1 ht 

«Xl(7,w) = + - - uj Xo 7 ), (36 

7 7 1 — 7 + u; 1 — 7 

where A\{oj) is related to the Mellin transform of the DGLAP splitting function and xjf i s the 
higher twist part of xo> 

Xo*(7) = Xo(7) - - - = 2^(1) - V(l + 7) " ^(2 - 7)- (37) 

7 1-7 
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0.2 0.4 0.6 0.8 1 



FIG. 2: Characteristic functions w(7, a = 0.2). The solid line is LL BFKL, Eq. ([2]), the dotted line is the 
asymmetrical kinematical constraint, Eq. (|34|) . and the dashed line is the KMRS model (|35l) . 

We make the approximation that np = in the evolution, which gives A±(u) ~ —11/12. It is 
possible to include the effects of quarks [l3| but they will be neglected here. 

The characteristic kernel (I35p is shown in Fig. [2] together with the LL BFKL kernel and the 
asymmetric consistency constraint. An important feature of the kernel (|35p is that it has no pole 
at 7 = 0. Instead it approximately fulfills energy conservation through 7(0; = 1) ~ (exact energy 
conservation implies equality). 



D. Comparison of results 

Let us now compare the results obtained from the different suggested kernels. We will focus on 
two quantities, the derivative of the saturation scale and the reduced front shape. 
The logarithmic derivative 

MY) - ^ (38) 

gives a measure of how fast the fronts propagate toward large momenta and should asymptotically 
approach the constant value v*. In Fig. O the results for X S (Y) obtained from the three different 
kernels are plotted together with the LL BFKL result. It is clear that the NLL corrections decrease 
v* , but also that the approach to the asymptotics is somewhat smoothed out. This can be traced 
to the smaller values of the diffusion constants. For experimentally observable rapidities it seems 
that the velocity is a decreasing function of Y. Remember, however, that for smaller rapidities 
there can be appreciable effects coming from the non-universal terms that are not included here 
(sec [5] for numerical investigations of this). 

The reduced front shape is obtained as the combination 

M(k, Y) = M(k/Q S (Y), Y) x (k 2 /Q 2 (Y))^ (39) 
1 ^(t^)^(-^ 2 (t£^)) (40) 



VD \Q 2 s (Y)J P V 4£>F \Q 2 s(Y) 
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FIG. 3: Derivative of the saturation scale, X S (Y), as defined in Eq. (|38|) for LL BFKL, the kinematical 
constraint (KC), the rapidity veto, and the KMRS kernel. 



M(k,Y) 
3 



KMRS 




ln(k 2 /Q s 2 ) 



FIG. 4: Reduced front shape as defined in Eq. (|4"0)) . obtained from LL BFKL, the kinematical constraint, 
and the KMRS kernel, and plotted for rapidities Y = 2,4, 6, 8, 10. 



which removes the exponential part and quantifies the deviation from geometric scaling. For perfect 
geometric scaling the reduced front should not depend on the rapidity. In Fig. [4] the reduced fronts 
are plotted for the LL BFKL kernel, the consistency constraint and the KMRS model. The curves 
have widths which reflect the diffusion constants in the diffusive breaking of geometric scaling, and 
it is obvious that the RG-improved kernels lead to a slower diffusive breaking of geometric scaling. 



VI. DIFFUSIVE DIFFERENTIAL EQUATION APPROXIMATION 



In l_9f] the BK equation was approximated as a partial differential equation, which could be 
transformed into the FKPP equation by expanding the BFKL kernel to second order. The RG- 
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improved BK equation may similarly be approximated as a differential equation, which will be a 
generalization of the FKPP equation. It is instructive to study this approximation since we then 
obtain closed analytical expressions and the structure of the modifications to the traveling wave 
equations may be easier appreciated. 

The BFKL function \ can be Taylor expanded around the point (7 C , tu c ) to second order, 

x(7> w ) ~ X + x\l - 7c) + x(w - w c ) 

+ \ [x"(7 " 7c) 2 + " ^c) 2 + 2 X '(7 " 7c) (w - w c )] (41) 

where in this Section the subscript c on x is omitted for simplicity of notation, so that \ = 
X(7c,w c ), x' = d~/X(lc,uJc), X = &>x(7c,w c ), x' = 9 7 aa;x(7c,w c ), and so on. The BK equation for 
J\f = J\f(L, Y) in the diffusive approximation then takes the form 



dyM 



a 



Y d L + ^ d Y + (X'lc ~X' + X'^c) d L - (x^ c -X + X7c) dy 



./ ., ., / / x"lo , ./ , X^c 2 



-X d L d Y + X - X 7c - X^ c H h X 7c^ c + 



N - aN 2 (42) 



which at first glance looks more symmetric in L and Y than the LL approximation. The equation 
contains a cross term derivative and a term proportional to 8l, but the variable change 

( i = AF + Bi (43) 
[x = CY + DL 

gives an equation on the form 

r dfu + d t u = d 2 x u + u - u 2 . (44) 

The explicit expressions for the constants r, A, B, C, and D in terms of x>x'>X e ^ c - are n °t very 
illuminating and are given in the Appendix. The function u is defined as 

/ x ,r( T Ax-Ct ^ Dt-Bx\ 

«'< t] = const x ( ArrrBc' Y = UTTbc J ■ (45) 

Remember that all the coefficients are just pure numbers once the saddle point values 7 C and lu c 
have been found. 

Eq. (I44p is an extended FKPP equation with a higher time derivative. Note that in LL BFKL the 
evolution occurs in rapidity, since the time variable t cx Y, but in the present case the evolution 
instead occurs in t = AY + BL. The w-dependent kernel thus mixes evolution in rapidity and 
momentum. 

We can now analyze Eq. (|44|) in the same way as we analyzed the full linearized BK equation 
in the LL or NLL cases. The result is 

v* = . 2 (46) 
\* = ^TT4^ (47) 

D - (T+iTr < 48 > 

This is to be compared with the results for the first-order in time FKPP equation, v* = 2, A* = 
1,D = 1. Taking into account the transformation back to the variables L and Y we again see that 
the asymptotic velocity is reduced compared to the LL value. 
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It turns out, however, to be much simpler to accept a "convection term" proportional to d x in 
the equation. This makes it possible to keep the t oc Y property and the change of variables is 
much simpler since the only constraint is to get rid of the mixed derivative term instead of getting 
rid of several terms. In fact, such a convection term will only change the velocity by an additive 
constant. 

Thus, if we instead in Eq. (|42|) make the variable change 

\ t = AY (49) 
y x = CY + DL, 

we get the transformed equation 

r d\u + dtu = d 2 j,u + a d x u + u — u 2 . (50) 

with 

/ x Ct t\ 
u(x,t) = const x AM L = — - — ,Y = -J . (51) 

The coefficients r, a, A, C, and D are again given in the Appendix. This equation has an additional 
linear derivative in x, with a "convection coefficient" a. With this change of variables we therefore 
see the problem in a different way; the evolution in time is now evolution in rapidity as in LL 
instead of in rapidity and momentum as with the earlier variable change. There is a higher order 
time derivative slowing down the spreading, as well as a convection term damping the wave. 
Unfortunately the expressions for v*, A* and D are now quite complicated and I choose to not 
write them here. 

The interesting points here are that when passing to RG-improved kernels, the spreading velocity 
decreases, and the equation becomes more symmetric in L and Y. In terms of linearized partial 
differential equation approximations, the LL BFKL equation is first order in time and second order 
in space while the DGLAP equation is first order in space and second order in time. The RG- 
improved equation which contains higher order corrections is instead second order in both time and 



space. This duality is explicit in Eqs. flUD and (|50|) . reminding one of the approach of Ref. [10(] 



VII. SUMMARY AND OUTLOOK 

I have proposed a simple way to use the traveling wave framework to compute the energy 
dependence of the saturation scale for RG-improved BFKL kernels with fixed QCD coupling. This 



completes the similar computation done for running coupling in [24]]. Higher order corrections 
are found to suppress both the increase with energy of the saturation scale and the speed of the 
diffusive spread of the front that breaks geometric scaling. The results obtained from three of the 
different kernels proposed in the literature are somewhat different in size. I have also studied the 
differential equation approximation of the evolution equation and found that the RG-resummation 
qualitatively changes the obtained equations. 

The study of RG-improvements to the BK equation at fixed or running coupling is the first 
step in understanding higher order corrections to saturation equations. Recent studies of how 
to rigorously include a running coupling into the BK and Balitsky-JIMWLK equations are not 



conclusive [141 ]. and the full NLL corrections to the BK equations have not been computed. It is 
not clear if the linearized version of the NLL BK equation would be the same as the NLL BFKL 
equation, but it is a reasonable assumption. Therefore, at this time, it is useful to consider the 
effect of RG-resummed kernels on the solutions of the BK equation. 
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It would be interesting to compare the results obtained in this paper to numerical simulations 
for various kernels. The results could also be used in comparing to experimental data along the 



lines of e.g. 39f] . Another interesting question is to ask what happens when fluctuations |26j are 



included in the picture. 
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APPENDIX: VARIABLE CHANGES FOR DIFFUSIVE APPROXIMATION 

The first change of variables, Eq. (143 jl . is defined by 



t = AY + BL 
x = CY + DL 



where 



(A.l) 



A = ax' 2 u c + xax" ~ x" ~ ax'x' ~ ax"x^c B > A ^ 

ax' 2 lc ~ axx' + X 1 + «X'X ~ ax"X7c 

D = X-x'lc-X^-l/a ^ (A3) 

X ~ X 7c - X U>c 

and where B and C are given by B = P4/P2 and C 2 = fy/fy with 

Pl = _l X ix' a (V ~ X' u c ) + x" (~X a + xu c a + I)) 2 | 
2 \ ("7cX' 2 - Xax' + X' + Xa(x' ~ x" lc)f 
y,, + 2x' {x'ajx'uc ~ x') +X" (xa~ Xu c a- 1)) \ 
a 7c X' 2 ~ X a x' + X' + X " (x' ~ x" 7c) / 



n I 112 ■■ -2 2 , // / ■ 2 -2 -12 2-2,-2 2-2 

P2 = [X X a 7c + X (X a - x 7c a + % ^ c a - 

2 X X 7c « 2 + 2 x X 7c ^c a 2 + 2 x lo c a - 2 x (x a w c + 1) a + 1) + 

a (x' - x' w c ) (2 a 7 C x' 2 + (-2 x " + X w c a + 2) x + XX «)) 

/ (a (a 7c x' 2 - X « X' + X' + X « (x' - x" 7c))) (A.5) 

a _ l_ ( X" '(~X« + X' '7ca + X^cQ + I) 2 
2 ^ ^(-X' + X'^c + X'^c) 2 
2x / (-X« + X / 7c« + X^c«+ 1) 
a (x' - X"7c - X'u c ) 

fa = i(x"7 2 -2x , 7c + 2x / Wc7c + XWc 2 + 2x-2x^c). (A.7) 

Finally, r is given by r = Pi 04/(3%. 

The second variable change, Eq. (|49p . is defined by the somewhat simpler form 

it = AY 
\x = CY + DL, 



+ X (A.6) 



(A. 



where 
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A=f, C = d€, D 2 = ^, (A.9) 

P2 X P3 



and where the coefficients /3j are 

A--f 



#2 = X^c + X 7c - X + 1/a 

P4 = - [X Xlc - X Ic + XX - x x - — 
XV " 

a X 7c i •/ i X / 

A> = — ^ h X 7c + X 7c - X w c + X- 



We also have r = f3\j3^/ ' (5\ and cr = yj fi^Jfhfe. 



[1] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rept. 100, 1 (1983). 

[2] A. H. Mueller, in proc. QCD perspectives on hot and dense matter, Eds. J. -P. Blaizot and E. Iancu. 
(Kluwer, Dordrecht 2002), pp 45-72, |arXiv:hep-ph/0111244| ; 

E. Iancu and R. Venugopalan, in Quark Gluon Plasma 3, Eds. R.C. Hwa and X.N. Wang (World 

Scientific, Singapore 2004), pp 249, [arXiv:hep-ph/0303204] . 
[3] Y. V. Kovchegov, Phys. Rev. D 60, 034008 (1999). 
[4] I. Balitsky, Nucl. Phys. B 463, 99 (1996); 

Phys. Rev. Lett. 81, 2024 (1998); 

Phys. Rev. D 60, 014020 (1999); 

Phys. Lett. B 518, 235 (2001). 
[5] R. Enberg, K. Golec-Biernat and S. Munier, Phys. Rev. D 72, 074021 (2005). 
[6] L. N. Lipatov, Sov. J. Nucl. Phys. 23, 338 (1976); 

E. A. Kuraev, L. N. Lipatov and V. S. Fadin, Sov. Phys. JETP 44, 443 (1976); 

45, 199 (1977); 

I. I. Balitsky and L. N. Lipatov, Sov. J. Nucl. Phys. 28, 822 (1978). 
[7] V. S. Fadin and L. N. Lipatov, Phys. Lett. B 429, 127 (1998); 
M. Ciafaloni and G. Camici, Phys. Lett. B 430, 349 (1998); 

G. Camici and M. Ciafaloni, Phys. Lett. B 412, 396 (1997) [Erratum-ibid. B 417, 390(1998)]. 
[8] G. P. Salam, Acta Phys. Polon. B 30, 3679 (1999). 
[9] G. P. Salam, J. High Energy Phys. 9807, 019 (1998); 

M. Ciafaloni and D. Colferai, Phys. Lett. B 452, 372 (1999); 

M. Ciafaloni, D. Colferai and G. P. Salam, Phys. Rev. D 60, 114036 (1999). 
[10] R. D. Ball and S. Forte, Phys. Lett. B 405, 317 (1997); 

G. Altarelli, R. D. Ball and S. Forte, Nucl. Phys. B 575, 313 (2000); 

Nucl. Phys. B 599, 383 (2001); 

Nucl. Phys. B 621, 359 (2002); 

Nucl. Phys. B 674, 459 (2003). 
[11] J. Kwiecihski, A. D. Martin and P. J. Sutton, Z. Phys. C 71, 585 (1996); 

B. Andersson, G. Gustafson, H. Kharraziha and J. Samuelsson, Z. Phys. C 71, 613 (1996); 

B. Andersson, G. Gustafson and J. Samuelsson, Nucl. Phys. B 467, 443 (1996). 
[12] C. R. Schmidt, Phys. Rev. D 60, 074003 (1999); 

J. R. Forshaw, D. A. Ross and A. Sabio Vera, Phys. Lett. B 455, 273 (1999). 
[13] V. A. Khoze, A. D. Martin, M. G. Ryskin and W. J. Stirling, Phys. Rev. D 70, 074013 (2004). 



17 



[14] E. Gardi, J. Kuokkanen, K. Rummukainen and H. Weigert, arXiv:hep-ph/0609087; 

Y. V. Kovchegov and H. Weigert, arXiv:hep-ph/0609090; 

I. Balitsky, |arXiv:hep-ph/0609105[ ~ 
[15] A. H. Mueller, Nucl. Phys. B 558, 285 (1999). 

[16] E. Iancu, K. Itakura and L. McLerran, Nucl. Phys. A 708, 327 (2002). 

[17] A.H. Mueller and D.N. Triantafyllopoulos, Nucl. Phys. B 640, (2002) 331. 

[18] D.N. Triantafyllopoulos, Nucl. Phys. B 648, 293 (2003). 

[19] S. Munier and R. Peschanski, Phys. Rev. Lett. 91, 232001 (2003). 

[20] S. Munier and R. Peschanski, Phys. Rev. D 69, 034008 (2004). 

[21] S. Munier and R. Peschanski, Phys. Rev. D 70, 077503 (2004). 

[22] A. M. Stasto, K. Golec-Biernat and J. Kwiecihski, Phys. Rev. Lett. 86, 596 (2001). 

[23] R. Enberg, talk given at the PIS 2005 workshop, Ap ril, 2005. 
[24] R. Peschanski and S. Sapeta, [arXiv:hep-ph/0610354[ 
[25] Y. V. Kovchegov, Phys. Rev. D 61, 074018 (2000). 

[26] E. Iancu, A. H. Mueller and S. Munier, Phys. Lett. B 606, 342 (2005). 

[27] R.A. Fisher, Annals Eugen. 7, 355 (1937), available at http:/ /www. library.adelaide.edu.au/digitised/fisher/152.pdf 
[28] A.N. Kolmogorov, I.G. Petrovskii, N.S. Piskunov, Bull. Univ. Etat a Moscou, Ser. int., A 1, 1 (1937), 

reprinted in English in P. Pelce, Dynamics of Curved Fronts (Academic Press, San Diego, 1988). 
[29] For recent reviews, see e.g. 

W. van Saarloos, Phys. Rept. 386, 29 (2003); 

J.D. Murray, Mathematical Biology, 3rd Ed., Springer Verlag (2002). 
[30] U. Ebert, W. van Saarloos, Physica D 146, 1 (2000). 
[31] M. Bramson, Mem. Am. Math. Soc. 44, No. 285 (1983). 

[32] K. Golec-Biernat, L. Motyka and A. M. Stasto, Phys. Rev. D 65, 074037 (2002). 
[33] G. Chachamis, M. Lublinsky and A. Sabio Vera, Nucl. Phys. A 748, 649 (2005). 
[34] E. Gotsman, E. Levin, U. Maor and E. Naftali, Nucl. Phys. A 750, 391 (2005). 
[35] K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 59, 014017 (1999); 

K. Golec-Biernat and M. Wiisthoff, Phys. Rev. D 60, 114023 (1999). 
[36] E. Brunet and B. Derrida, Phys. Rev. E 56, 2597 (1997). 

[37] R. Enberg, B. Pire, L. Szymanowski and S. Wallon, Eur. Phys. J. C 45, 759 (2006). 

[38] D. Y. Ivanov and A. Papa, Nucl. Phys. B 732, 183 (2006). 

[39] E. Iancu, K. Itakura and S. Munier, Phys. Lett. B 590, 199 (2004). 



